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Applying the phenomenon of neutrino lasing in the solar interior, we 
show how the rate for the generic neutrino decay process u — fermion + 
boson, can in principal be enhanced by many orders of magnitude over 
its normal decay rate. Such a large enhancement could be of import to 
neutrino- decay models invoked in response to the apparent deficit of elec- 
tron neutrinos observed from the sun. The significance of this result to 
such models depends on the specific form of the neutrino decay, and the 
particle model within which it is embedded. 



PACS number(s):12.15.Cc, 13.35.+S, 14.60. Gh, 96.60.Kx 



1 



The Solar Neutrino Problem (SNP) is well established. Neutrino fluxes predicted 
by standard solar models [1] are in disagreement with the fluxes observed by current 
experiments. In terms of the ratio of observed-to-predicted neutrino fliixes, rj, the Cl- 
Ar experiment [2] reports rjci = 0.28 ± 0.03; whereas the Kamiokande collaboration [3] 
report ijic = 0.51 ± 0.07. These two experiments are mostly sensitive to the high-energy 

neutrinos. In order to measure the lower-energy pp neutrinos, two Ga-Ge experiments 
have been attempted. The SAGE collaboration [4] report rjsAGE = 0.53 ± 0.19; whereas, 
the GALLEX collaboration [5] report ijgall = 0.62±0.15. The errors are the la statistical 
and systematic errors reported by the diff'erent groups added in quadrature: theoretical 
uncertainties from the standard solar model are not included. 

Neutrino decay is often invoked as a possible remedy to the the above discrepancies 
[6,7]. Although Ue decay as a solution to the SNP would appear to be in contradiction with 
the observed z>e pulse from SN1987A, a viable solution remains where the mass eigenstates 
vi ,1/2 are a substantial mixture of the weak eigenstates Ve ,i^x {x = ^,r) [8]. Matter 
induced effects can also require some additional interpretation of the supernova limits [9] . 
The generic neutrino decay can be described V2 z^i + -S, where B is some unknown 
boson. The most favored identification of the boson is with some type of majoron particle 
(f). For example, the two decay modes V2 ^ vi + (f) and V2 ^ + have been previously 
analysed in detail, and regions of the lifetime and mixing-angle phase-space consistent with 
current observations determined [7]. Fast majoron decay into sterile neutrino states has 
also been discussed with regard to the SNP [10]. 

The main purpose of this report is to point out that the phenomenon of neutrino lasing 
may be of import with regard to the generic neutrino decay paradigm. Specifically, we 
will find that the bosons emitted in relativistic neutrino decays can stimulate the decay to 
proceed orders of magnitude faster than that described by the normal decay rate. Regions 
in the allowed parameter space consistent with the observed neutrino fiuxes are therefore 
subject to revision. 
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The phenomenon of neutrino lasing in the context of the early universe and the dark- 
matter problem has previously been discussed [11]. Neutrino lasing can best be described 
as a process in which the decay of a relativistic neutrino proceeds by stimulated emission of 
bosons, thereby dramatically increasing its decay rate. We wish to see how solar neutrinos 
can have their decay rates affected by this phenomenon. Considering the generic decay 
(and its inverse) of a heavy neutrino, i/, into a lighter neutrino, F and a boson S, it can 
be shown that [11] the evolution of the occupation number distributions /s, is described 
by the following Boltzmann equation: 

fB{EB) = g ° / dEHiMl - IfKI + Ib) - /b/f(1 - Ih)] (1) 

moLBPE Je-{Eb) 

where is the boson three- momentum, Tq is the free decay rate for a at rest, and 
mo/2 is the three- momentum of the decay products in the H rest- frame. Here we have 
assumed isotropic decay in the H rest frame. In the massless boson limit we have 



mo = mn ■ (2) 



The integration is over the energy-conserving plane E}j — Ep + Eb, with limiting values 
and 

Eb,f = imo/2mH) {Eh ± Ph) ■ (36) 
Considering first only the initial growth of /s, our Boltzmann equation becomes 



h ~ ^/sro / dEfn , (4) 



where we have taken the growth rate to be maximized at (^i mjj/EH)- We can 
approximate 

dEfn - . (5) 



Therefore, 

^ -Eff 3 
rriH 

For significant decay of particles we require p%fB ~ PnfH, that is we need 



fB - (^) . (7) 



rrijj 



Therefore from eqs (6) and (7) we find that the time required for the number of bosons to 
equal the number of neutrinos at some point in the sun is given by 



^f^)'ln(^)l, (8) 



which typically is much less than the age of the sun. 

Using Fo = {Eh/uih)^®, where F© is the decay rate in the solar frame, Eq. (6) allows 
us to define a new efi^ective decay rate in the presence of lasing as 



lase 



^f^)>o. (9) 



As an example of the effect consider the decay of the pp neutrinos as they pass through 
the point R = O.li?©. Prom the pp neutrino density at O.li?©, we find ^ ~ 10~^^, and 
therefore for Ejj/mH ^ 10'^, the decay rate is enhanced. For example, if the mass of 
the electron neutrino was uih ~ 0.01 eV, the lasing rate is a factor ~ 10^ larger than 
F0. Clearly for even smaller values of the neutrino mass the effect on the decay rate will 
become more significant. This would be the case as long as E^ ^ since if this bound 

4 



is not satisfied the above analysis would break down as the wavelength of the emitted 
boson would be larger than the size of the region producing it. For practical purposes, 
we can impose the limit mn > lO"^** eV as the mass above which our analysis would be 
valid. Another limitation of our analysis is that we have neglected transport terms in the 
Boltzmann equation. This is a good approximation as long as Tiase » since then 
fi) will be able to grow. Finally, note from the / terms of Eq. (1) that the lasing effect 
"switches off" when the abundance of decaying neutrinos equals that of the decay-product 
neutrinos. 

As mentioned earlier, we have assumed in the above analysis that decays in the rest 
frame of the decaying neutrino are isotropic. What form does Eq. (9) take if we drop this 
assumption? In the more general case, the Boltzmann equation Eq. (1) takes the form [12] 



fB = g ° / dEHlMl - /f)(1 + Ib) - - m 

moLBPB JE-{Eb) 



X 



1 ^^E^ + E-^-'^Eb 



rriH PHrrio 



(10) 



where a is defined through the probability distribution P{9) given by 



P(^) = l(l-acos^), (11) 



and where 9 is the angle in the rest frame between the decay product velocity and the 
parent velocity in the solar frame. The isotropic case discussed previously corresponds to 
q; = 0, and is valid (up to factors of ~ 2) for all a > —1. For a ~ — 1, the emission of 
low-momentum bosons is suppressed. Utilizing Eq. (10), and applying similar arguments 
as those used in the isotropic case we find the equivalent equation to Eq. (9) is 

-L lase ~ T^Q- J J- © 
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(12) 



From comparison of Eq. (9) with (12) it can be seen that a = —1 results in a suppression 
factor {mH/Eff)'^, relative to the case of isotropic decays. 

Let us consider the pp neutrinos produced by the sun, and focus on the decay of these 
neutrinos into a tau or muon neutrino {vg Vx + B) . In order to accurately assess the 
importance of neutrino lasing we must average the effect over the varying matter and 
neutrino densities in the sun. Let us assume that matter effects dominate the neutrino 
masses. In this case we can replace uih in Eqs. (8) and (12) with y/2EV, where E is the 
Ue energy and the effective potential V is given by 

V = V2GFpm-^Ye . (13) 

Here, Gp is the Fermi constant, p is the matter density, rrin is the nucleon mass, and Fg is 
the number of electrons per nucleon. 

In the solar frame, the decay rate for some process described by some dimensionless 
coupling constant g can be written 



^® - 16^ = 16^ • 



Neglecting lasing effects, in terms of the electron-neutrino flux at some origin, the flux 
^e{r) at distance r from the origin is given by [13] 



$e(r) = K exp E 9^^""] ' (15) 



where the effective matter width d^'^ ^ for the decay is given by 



d^ff = / p{r)Ye{r)dr, (16) 







and where the refraction width do is given by do = \/2'kG p^nin- 

The reduction in the neutrino flux when lasing effects are included can be approximated 

by 

^e{rf = K^^-^Y.9^D^ff] , (17) 



8do ^ 



where 



and where 



i:)^// ~ f x{r)p{r)Ye{r)dr , (18) 
Jo 



71 

= EVW) ■ 

We apply Eq. (17) to a solar model [14], and determine the ratio ^e{f)/^e{fy°'^^- 
Our calculations show that for a decay such as ^-g — ^^a; + where the chirality of the 
neutrinos remain unchanged, $e('')/^e('')'"*^ can be >> 1 for both pp and neutrinos. 
This means that any such decay applied to the solar neutrinos must take into account this 
effect. For processes where a = —1, for example chirality flipping decay with spin-zero 
boson emission, Eq. (19) becomes 



and there is a negligible effect due to the spin suppression effects discussed earlier. 

In order to see the effects of lasing let us assume a chirality-conserving interaction 
and apply the current laboratory bounds on the chirality-conserving decay t-'x + 4>^ 

namely g as defined in Eq. 14 is < 7.0 x 10""^ [15]. Remembering that the lasing switches 
off when fy^ = fy^ , we find that half the pp neutrinos could in fact undergo decay in the 
solar interior. This conclusion remains valid for g ^ 10~^. Decays of Vf. to a left-handed 
antineutrino through majoron emission give a similar result (though a slight modification 
to the potential V is required in this case). Decay into a right-handed antineutrino would 
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require that the chirahty be flipped, and thus is unimportant due to spin suppression 
effects. 

A question remains as to whether there is a weU-defined model which couples neutrinos 
to a boson in such a way that we will see lasing. As stated above the most obvious candidate 
to consider is the singlet majoron model [16]. Consider a scalar potential of the form 

V[^] = A($2 _ ^2)2 _ (21) 

Clearly the vacuum corresponds to |($)| = v, but we have a number of ways of describing 
excitations. If we consider the linear realization of the symmetry breaking, then we write 

^ = v + p + i(f) (22) 

and find that p is a massive scalar field, and that is a massless Goldstone boson - the 
majoron. Including matter effects, we can determine that part of the lagrangian which 
involves the neutrinos and the linear majoron coupling: 



NrMdNl + \NRMR{NRy + iy.c] + ^^NRMRiNRf - I^C^NlI^Nl , (23) 



where N^r^ are the weak-lepton-doublet (singlet) neutrinos, (^u = GFP'm~^YnVn (where 
Vfj_ is collective nucleon four- velocity), and the mass matrices Md and Mr are assumed 
real. Considering only two generations and assuming Mr » Mjj, diagonalization of the 
mass matrix determines the physical fields ui^2- The coupling of these neutrino fields to 
in this case is direct, and has the form 
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(24) 



meaning that the local operator for the transition 1/2 — > i'\-\-4> necessarily involves a chirality 
flip. As seen earlier, the spin suppression associated with this chirality flip renders lasing in 
the sun unimportant. Only within environments where the neutrino density is much larger 
than solar neutrino densities (eg. supernova explosions), would the coupling of Eq. (24) 
be important. 

It is possible to have chirality-conserving decays from other processes. In the linear 
expansion this occurs when mass insertions are added to the external neutrino lines. How- 
ever, in the most simple models each mass insertion introduces a suppression factor of 
order (m/i?)^. If the exponential expansion $ = (f + p) exp{i(j) / v) is used, the chirality 
conserving decays naively involve derivative couplings which suppress the emission of low 
momentum bosons, and hence the lasing rate also. The actual result must be independent 
of which expansion is utilized, but only after all relevant diagrams at each order are in- 
cluded would the exact suppression factor be determined. Such suppression would clearly 
be model dependent. Indeed, we note that it is possible to construct more complicated 
majoron models which exhibit no suppression of the chirality-conserving decays [18]. 

Another hopeful possibility is that the neutrino decays through a new spin-1 field. The 
coupling would be between neutrinos of the same chirality state, and the rate would not 
be suppressed by any chirality or momentum effects. Massive vector fields were speculated 
on as possible sources of lepton number violation through local symmetry breaking [17], 
but were superseded by majoron models and global symmetry breaking, which seemed to 
arise more naturally in GUT's. 

In summary, we have seen how neutrino lasing may play a role in neutrino decay 
processes occurring in the solar interior. We have highlighted how the spin nature of the 
particles involved in the decay process play an important role. In particular, we noted that 
within the context of the most simplest singlet majoron models, no significant neutrino 
lasing in the sun would proceed. More complicated extensions to the standard model need 
to be invoked if lasing is to be viable in the solar interior. Since the energy spectrum of the 
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decaying electron neutrinos is not degraded by the lasing phenomenon, future detectors 
should be able to distinguish between the process described here and that anticipated from 
the usual decay paradigm. 

A similar analysis of neutrino lasing applied to supernovae explosions may provide 
some additional constraints. We note again, however, that since bosons are preferentially 
produced at low momentum the energetics of such explosions should not be significantly 
altered by excess energy losses, unless the decay product neutrino is sterile. 

Neutrino lasing has the potential to be an important phenomenon. In addition to its 
impact on the dark matter problem, it may also have interesting implications for majoron 
models of baryogenesis due to enhanced decays of sterile neutrinos. Any experimental 
information on the viability of neutrino lasing would therefore be very valuable. 

We thank N. Kaiser, M. Luke, S. Selipsky and L. Widrow for useful discussions. 
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